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2  Publications 


The  following  publications  all  acknowledge  ARO  Grant  DAAH04-96-1-0013.  For  brevity, 
we  do  not  list  papers  published  in  conference  proceedings  for  which  a  journal  version  is  also 
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[1]  P.  K.  Agarwal.  Range  searching.  In  J.  E.  Goodman  and  J.  O’Rourke,  editors,  Handbook 
of  Discrete  and  Computational  Geometry,  chapter  31,  pages  575-598.  CRC  Press  LLC, 
Boca  Raton,  FL,  1997. 

[2]  P.  K.  Agarwal,  N.  Amenta,  B.  Aronov,  and  M.  Sharir.  Largest  placements  and  motion 
planning  of  a  convex  polygon.  In  J.-P.  Laumond  and  M.  Overmars,  editors,  Algorithms 
for  Robotic  Motion  and  Manipulation,  pages  143-154,  Wellesley,  Massachusetts,  1997. 
A.K.  Peters.  Proc.  1996  Workshop  on  the  Algorithmic  Foundations  of  Robotics. 

[3]  P.  K.  Agarwal,  N.  Amenta,  and  M.  Sharir.  Placement  of  one  convex  polygon  inside 
another.  Discrete  Comput.  Geom.,  19:95-104,  1998. 

[4]  P.  K.  Agarwal,  L.  Arge,  G.  Brodal,  and  J.  S.  Vitter.  1/ o-efficient  dynamic  point  location 
in  monotone  subdivisions.  In  Proc.  10th  ACM- SIAM  Sympos.  Discrete  Algorithms, 
pages  11-20,  1999. 

[5]  P.  K.  Agarwal,  L.  Arge,  and  J.  Erickson.  Indexing  moving  points.  In  Proc.  Annu. 
ACM  Sympos.  Principles  Database  Syst.,  2000.  175-186. 

[6]  P.  K.  Agarwal,  L.  Arge,  J.  Erickson,  P.  G.  Franciosa,  and  J.  S.  Vitter.  Efficient 
searching  with  linear  constraints.  Journal  of  Computer  and  System  Sciences,  pages 
194-216,  October  2000. 

[7]  P.  K.  Agarwal,  L.  Arge,  T.  Murali,  K.  Varadarajan,  and  J.  Vitter.  I/O-efficient  al¬ 
gorithms  for  contour  line  extraction  and  planar  graph  blocking.  In  Proc.  ACM-SIAM 
Symp.  on  Discrete  Algorithms ,  pages  117-126,  1998. 

[8]  P.  K.  Agarwal,  L.  Arge,  O.  Procopiuc,  and  J.  S.  Vitter.  A  framework  for  index  bulk 
loading  and  dynamization.  In  Proc.  28th  Inti.  Conference  on  Automata  and  Program¬ 
ming  Langs.,  2001. 

[9]  P.  K.  Agarwal,  L.  Arge,  and  J.  Vahrenhold.  Time  responsive  indexing  schemes  for 
moving  points.  In  Proc.  7th  Workshop  on  Data  Structures,  Lecture  Notes  Comput. 
Sci.  Springer- Verlag,  2001. 

[10]  P.  K.  Agarwal,  B.  Aronov,  T.  M.  Chan,  and  M.  Sharir.  On  levels  in  arrangements  of 
lines,  segments,  planes,  and  triangles.  Discrete  Comput.  Geom.,  19:315-331,  1998. 

[11]  P.  K.  Agarwal,  B.  Aronov,  S.  Har-Peled,  and  M.  Sharir.  Approximation  and  exact 
algorithms  for  minimum- width  annuli  and  shells.  In  Proc.  15th  Annu.  ACM  Sympos. 
Comput.  Geom.,  pages  380-389,  1999. 
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[12]  P.  K.  Agarwal,  B.  Aronov,  J.  O’Rourke,  and  C.  A.  Schevon.  Star  unfolding  of  a 
polytope  with  applications.  SIAM  J.  Comput.,  26:1689-1713,  1997. 

[13]  P.  K.  Agarwal,  B.  Aronov,  J.  Pach,  R.  Pollack,  and  M.  Sharir.  Quasi-planar  graphs 
have  a  linear  number  of  edges.  Combinatorica,  17:1-9,  1997. 

[14]  P.  K.  Agarwal,  B.  Aronov,  and  M.  Sharir.  Line  traversals  of  balls  and  smallest  enclosing 
cylinders  in  three  dimensions.  In  Proc.  8th  ACM-SIAM  Sympos.  Discrete  Algorithms , 
pages  483-492,  1997. 

[15]  P.  K.  Agarwal,  B.  Aronov,  and  M.  Sharir.  On  levels  in  arrangements  of  lines,  segments, 
planes  and  triangles.  In  Proc.  ACM  Symp.  on  Computational  Geometry ,  pages  30-38, 
1997. 

[16]  P.  K.  Agarwal,  B.  Aronov,  and  M.  Sharir.  Motion  planning  for  a  convex  polygon  in  a 
polygonal  environment.  Discrete  Comput.  Geom.,  22:201-221,  1999. 

[17]  P.  K.  Agarwal,  B.  Aronov,  and  M.  Sharir.  Exact  and  approximation  algorithms  for 
minimum-width  cylindrical  shells.  In  Proc.  11th  ACM-SIAM  Sympos.  Discrete  Algo¬ 
rithms,  pages  510-517,  2000. 

[18]  P.  K.  Agarwal,  B.  Aronov,  and  M.  Sharir.  On  the  complexity  of  many  faces  in  arrange¬ 
ments  of  circles.  In  Proc.  41st  Annual  IEEE  Sympos.  on  Foundations  of  Computer 
Science,  2001. 

[19]  P.  K.  Agarwal,  J.  Basch,  L.  Guibas,  and  J.  Hershberger.  Lower  bounds  for  kinetic 
planar  subdivisions.  Discrete  Comput  Geom.,  24:721-733,  2000. 

[20]  P.  K.  Agarwal,  J.  Basch,  L.  J.  Guibas,  J.  Hershberger,  and  L.  Zhang.  Deformable  free 
space  tiling  for  kinetic  collision  detection.  In  Proc.  4th  Workshop  Algorithmic  Found. 
Robot.,  2000.  To  appear. 

[21]  P.  K.  Agarwal,  T.  Biedl,  S.  Lazard,  S.  Robbins,  S.  Suri,  and  S.  Whitesides.  Curvature- 
constrained  shortest  paths  in  a  convex  polygon.  In  Proc.  14th  Annu.  ACM  Sympos. 
Comput.  Geom.,  page  to  appear,  1998. 

[22]  P.  K.  Agarwal,  M.  de  Berg,  J.  Gudmundsson,  M.  Hammar,  and  H.  J.  Haverkort.  Box- 
trees  and  r-trees  with  near-optimal  query  time.  In  Proc.  17th  Annu.  Symposium  on 
Computational  Geometry,  2001. 

[23]  P.  K.  Agarwal,  M.  de  Berg,  D.  Halperin,  and  M.  Sharir.  Efficient  generation  of  k- 
directional  assembly  sequences.  In  Proc.  7th  ACM-SIAM  Sympos.  Discrete  Algorithms, 
pages  122-131,  1996. 

[24]  P.  K.  Agarwal,  M.  de  Berg,  S.  Har-Peled,  M.  Overmars,  M.  Sharir,  and  J.  Vahrenhold. 
Reporting  all  intersecting  pairs  of  polytopes  in  two  and  three  dimensions.  In  Proc.  7th 
Workshop  on  Data  Structures,  Lecture  Notes  Comput.  Sci.  Springer- Verlag,  2001. 
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[25]  P.  K.  Agarwal  and  P.  K.  Desikan.  An  efficient  algorithm  for  terrain  simplification.  In 
Proc.  8th  ACM-SIAM  Sympos.  Discrete  Algorithms ,  pages  139-147,  1997. 

[26]  P.  K.  Agarwal  and  P.  K.  Desikan.  Approximation  algorithms  for  layered  manufacturing. 
In  Proc.  11th  ACM-SIAM  Sympos.  Discrete  Algorithms,  pages  528-537,  2000. 

[27]  P.  K.  Agarwal,  D.  Eppstein,  L.  J.  Guibas,  and  M.  Henzinger.  Parametric  and  kinetic 
minimum  spanning  trees.  In  Proc.  39th  Annu.  IEEE  Sympos.  Found.  Comput.  Sci., 
pages  596-605,  1998. 

[28]  P.  K.  Agarwal  and  J.  Erickson.  Geometric  range  searching  and  its  relatives.  In 
B.  Chazelle,  J.  E.  Goodman,  and  R.  Pollack,  editors,  Advances  in  Discrete  and  Compu¬ 
tational  Geometry,  volume  223  of  Contemporary  Mathematics,  pages  1-56.  American 
Mathematical  Society,  Providence,  RI,  1999. 

[29]  P.  K.  Agarwal,  J.  Erickson,  and  L.  J.  Guibas.  Kinetic  BSPs  for  intersecting  segments 
and  disjoint  triangles.  In  Proc.  9th  ACM-SIAM  Sympos.  Discrete  Algorithms,  pages 
107-116,  1998. 

[30]  P.  K.  Agarwal,  E.  Flato,  and  D.  Halperin.  Polygon  decomposition  for  efficient  construc¬ 
tion  of  Minkowski  sums.  In  Proc.  8th  Annu.  European  Sympos.  Algorithms,  volume 
1879  of  Lecture  Notes  Comput.  Sci,  pages  20-31.  Springer- Verlag,  2000. 

[31]  P.  K.  Agarwal,  E.  F.  Grove,  T.  M.  Murali,  and  J.  S.  Vitter.  Binary  space  partitions 
for  fat  rectangles.  SIAM  J.  Comput.,  29:1422-1448,  2000. 

[32]  P.  K.  Agarwal,  L.  J.  Guibas,  S.  Har-Peled,  A.  Rabinovitch,  and  M.  Sharir.  Computing 
the  penetration  depth  of  two  convex  polytopes  in  3d.  In  Proc.  7th  Scandanavian 
Workshop  on  Algorithmic  Theory,  page  ??,  2000. 

[33]  P.  K.  Agarwal,  L.  J.  Guibas,  J.  Hershberger,  and  E.  Veach.  Maintaining  the  extent  of 
a  moving  point  set.  In  Proc.  5th  Workshop  Algorithms  Data  Struct.,  volume  1272  of 
Lecture  Notes  Comput.  Sci,  pages  31-44.  Springer- Verlag,  1997. 

[34]  P.  K.  Agarwal,  L.  J.  Guibas,  T.  M.  Murali,  and  J.  S.  Vitter.  Cylindrical  static  and 
kinetic  binary  space  partitions.  Computational  Geometry,  16(2):  103-127,  2000. 

[35]  P.  K.  Agarwal  and  S.  Har-Peled.  Maintaining  approximate  extent  measures  of  moving 
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2000. 
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3  Scientific  Progress  and  Accomplishments 

3.1  Introduction  and  Overview 

Geometric  computing,  which  emerged  as  a  self-standing  discipline  about  a  quarter  of  a  cen¬ 
tury  ago  under  the  broad  rubric  of  Computational  Geometry,  was  strongly  motivated  by 
the  pervasive  relevance  of  geometry  to  a  large  variety  of  applied  fields,  and  was  originally 
viewed  as  an  exciting  topic  in  the  design  and  analysis  of  algorithms.  In  the  context  of  this 
project,  we  underscore  that  Computational  Geometry  is  a  powerful  aid  to  a  number  of  mili¬ 
tary  and  civilian  applications  involving  a  physical  or  visual  environment,  including  robotics, 
navigation,  automatic  target  recognition  (ATR),  command,  control,  and  communications 
(C3),  battlefield  management,  geographic  information  systems,  computational  metrology, 
hit  avoidance,  computer-aided  design,  and  computer  graphics. 

Over  the  years  researchers  in  computational  geometry  developed  an  impressive  corpus  of 
algorithms  and  algorithmic  paradigms,  but  became  also  aware  of  the  severe  limitations  of  the 
original  model  based  on  the  assumed  availability  of  exact  real-number  arithmetic.  Instead, 
the  reality  of  computing,  which  results  in  nonrobust  implementations  of  mathematically 
correct  algorithms,  strongly  hampered  the  projected  transfer  of  technology  that  was  the 
fundamental  motivation  of  the  field. 

As  a  consequence  of  this  unsatisfactory  relation  to  real-world  applications,  in  the  last 
decade  the  field  underwent  a  significant  philosophical  transformation,  which  positioned  it 
strategically  for  the  materialization  of  the  desired  connection  with  the  applied  areas.  Specif¬ 
ically,  a  careful  study  of  robust  geometric  computing  and  the  serious  development  of  readily 
available  libraries  of  reliable  geometric  applications  were  the  two  most  significant  compo¬ 
nents  of  the  proposal  for  the  present  project,  which  was  submitted  in  response  to  the  ARO 
MURI-initiative  in  early  1995.  Funding  for  the  five-year  undertaking  was  granted  in  the 
summer  of  1995. 

At  the  completion  of  the  funded  research,  we  can  confidently  say  that  our  coordinated 
efforts  have  contributed  to  the  strengthening  of  the  discipline  and  to  the  forging  of  an  effective 
partnership  between  academic  research  and  practice,  both  civilian  and  military.  In  fact,  this 
re-orientation  of  the  research  community  [169]  has  been  met  by  a  complementary  acceptance 
by  the  application  community,  which  overcame  its  initial  skepticism  and  recognized  the 
invaluable  power  of  geometric  computing,  because  the  latter  is  grounded  on  solid  algorithmic 
and  topological  knowledge  and  uses  the  most  advanced  computing  technologies. 

In  summary,  we  believe  that  our  work  has  helped  prepare  a  most  fertile  ground  for 
the  further  evolution  of  the  field  in  a  rapidly  changing  technological  landscape.  Modern 
computers  have  processors  and  graphic  engines  that  are  fast  enough  to  perform  significant 
computations  on  collections  of  2-  and  3-dimensional  geometric  objects,  such  as  points,  line 
segments,  polygons,  curves,  surfaces,  etc.  In  addition,  the  advent  of  the  Internet  has  brought 
about  new  paradigms  and  modalities  for  using  computers  to  solve  large  problems,  which 
are  often  geometric.  Moreover,  these  paradigms  and  modalities  often  involve  interesting 
connections  to  other  areas  of  computer  science,  including  distributed  computing,  computer 
science  education,  and  electronic  commerce. 

Research  initiatives  and  detailed  individual  results  the  previous  five  interim  reports,  duly 
generated  at  the  conclusion  of  each  year  of  funded  research.  In  this  final  report,  however, 
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along  with  the  major  thematic  highlights  of  the  project,  we  emphasize  our  most  recent 
research. 

In  a  nutshell,  the  guiding  principles  of  our  research  have  been  the  pursuit  of  algorithmic 
robustness  and  ease-of-programming,  as  well  as  the  effective  management  of  scale.  These 
are  the  features  of  algorithm  engineering,  which  is  key  to  the  targeted  transfer  of  technology. 
Experimental  verification  was  an  important  component  of  this  effort  as  well. 

Our  presentation  is  organized  as  follows.  In  Section  3.2,  “Foundations”,  we  review  our 
contributions  to  basic  algorithmic  research  in  computational  geometry,  such  as  algorithms 
pertaining  to  geometric  structures  of  paramount  significance  (such  as  arrangements)  and 
general  data  organization.  In  Section  3.3,  “Methodologies”,  we  summarize  our  contributions 
to  correcting  the  model  inadequacies  of  traditional  computational  geometry,  specifically: 

•  methods  for  the  design  of  practical  and  robust  geometric  algorithms, 

•  development  of  a  finer  performance  framework  than  asymptotic  analysis, 

•  external  memory  management  techniques  for  large-scale  geometric  computations, 

•  a  transparent  parallel  I/O  environment  for  large-scale  geometric  computations  (TPIE), 

•  techniques  for  effectively  integrating  I/O  systems  into  massively  parallel  computers. 

In  Section  3.4,  “Emerging  Technologies”,  we  review  the  adaptation  of  our  research  effort  to 
the  new  computing  realities,  principally  new  powerful  software-engineering  paradigms  and 
the  Internet: 

•  an  object-oriented  library  of  robust  geometric  algorithms  (GeomLib), 

•  development  of  a  distributed  architecture  for  geometric  computing  over  the  Web  (Ge- 
omNet), 

•  development  of  approximation  algorithms  for  geometric  optimization  problems, 

•  algorithms  and  data  structures  for  mobile  data. 

Finally,  in  Section  3.5,  “Applications”,  we  discuss  our  numerous  contributions  to  a  variety 
of  applied  areas. 

3.2  Foundations 

Given  our  emphasis  on  understanding  the  fundamental  aspects  of  geometric  computing, 
we  studied  in  1999  several  fundamental  problems  in  combinatorics  and  in  computational 
geometry.  We  show  in  [76],  using  several  methods,  that  there  exists  a  set  S  of  n  points  in 
the  d-dimensional  unit  cube  so  that  every  d+ 1  points  of  S  define  a  simplex  of  volume  fi(^). 

In  [82]  we  define  a  new  type  of  distance  functions  in  the  plane  from  a  point  to  a  pair  of 
points.  We  focus  on  a  few  such  distance  functions,  analyze  the  structure  and  complexity  of 
the  corresponding  nearest-  and  furthest-neighbor  Voronoi  diagrams  (in  which  every  region 
is  defined  by  a  pair  of  point  sites),  and  show  how  to  compute  the  diagrams  efficiently. 
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A  graph  H  is  a  universal  graph  for  a  family  of  graphs  if  every  graph  in  the  family 
is  a  subgraph  of  H.  In  [107]  we  construct  a  universal  graph  of  size  0(n)  for  the  family 
of  n-vertex  planar  graphs  of  bounded  degree.  This  is  a  significant  improvement  over  the 
previously  known  bound  of  0(n  log n).  In  [108]  we  construct  a  universal  graph  of  size  0(n ) 
for  the  family  of  n-vertex,  bounded  degree  graphs  with  a  bisector  of  size  0(nc),  for  some 
c  <  1.  This  is  a  significant  improvement  over  the  previously  known  bound  of  0(n2c). 

3.2.1  Arrangements  of  surfaces 

We  studied  several  problems  involving  arrangements  of  surfaces.  In  [15,  10],  we  considered 
the  problem  of  bounding  the  complexity  of  the  A>th  level  in  an  arrangement  of  n  curves 
or  surfaces,  a  problem  dual  to,  and  an  extension  of,  the  well-known  k-set  problem.  Among 
other  results,  we  prove  a  new  bound,  0(nk5/3),  on  the  complexity  of  the  k- th  level  in  an 
arrangement  of  n  planes  in  Rz ,  or  on  the  number  of  Assets  in  a  set  of  n  points  in  three 
dimensions,  and  we  show  that  the  complexity  of  the  k-th  level  in  an  arrangement  of  n  line 
segments  in  the  plane  is  0(nVka(n/k)),  and  that  the  complexity  of  the  k- th  level  in  an 
arrangement  of  n  triangles  in  3-space  is  0(n2k^6a(n/k)). 

In  another  result,  we  proved  a  subcubic  bound  on  the  complexity  of  the  union  of  n 
congruent  cylinders.  This  is  the  first  step  toward  proving  a  near-optimal  bound  on  the 
complexity  of  the  Voronoi  diagram  of  lines  or  of  polygons  in  3-space. 

In  [18],  we  obtained  improved  bounds  on  the  complexity  of  many  faces  in  arrangements  of 
circles  and  unit  circles.  In  [52],  we  extended  the  duality  transform  to  a  family  of  pseudo- lines 
(a  family  of  curves,  any  two  of  which  intersect  in  at  most  one  point),  described  an  algorithm 
for  computing  this  duality,  and  showed  how  this  transform  can  be  used  to  develop  efficient 
algorithms  for  a  number  of  problems  involving  circles  and  other  quadratic  curves.  We  wrote 
two  comprehensive  surveys  on  this  topic  [49,  48]. 

3.2.2  Geometric  Algorithms  and  Data  Structures 

We  give  very  fast  work-efficient  parallel  algorithms  [138]  for  finding  convex  hulls  of  planar 
point  sets  and  solving  fixed-dimensional  linear  programs  [149]  in  a  computational  environ¬ 
ment  that  allows  for  fine-grain  parallelism.  Our  work  for  linear  programming  actually  gave  a 
general  method  for  efficiently  derandomizing  randomized  geometric  algorithms  sequentially 
or  in  parallel  and  uses  linear  programming  as  a  motivating  example.  We  also  give  optimal 
coarse-grain  parallel  algorithms  [141]  for  convex  hull  construction  and  a  general  class  of  data 
structure  searching  problems.  We  review  and  survey  [140]  in  an  invited  book  chapter  the 
major  results  in  parallel  computational  geometry,  many  of  which  are  the  results  of  Goodrich, 
Preparata,  Tamassia,  and  Vitter. 

We  define  and  study  a  combinatorial  problem  [156]  called  the  Weighted  Diagnostic  Cover 
(WDC),  which  models  genotyping  in  the  diagnosis  of  a  class  of  chromosomal  aberrations. 
We  develop  several  approximation  algorithms  for  WDC.  We  establish  worst-case  performance 
bounds  for  these  algorithms.  We  also  report  their  performance  on  a  real  data  set. 

In  [24],  we  developed  output-sensitive  algorithms  for  reporting  all  pairs  of  intersecting 
polytopes  in  a  given  set  of  convex  polytopes  in  3-space.  This  is  the  first  truly  output-sensitive 
algorithm  for  the  problem. 
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Motivated  by  the  fact  that  geometric  computing  naturally  involves  large  volumes  of 
pictorial  images,  and  considering  the  high  demand  of  raw  pictures  both  in  transmission  over 
the  net  and  in  local  storage,  we  have  thought  it  appropriate  to  revisit  the  effectiveness  of 
competing  data  compression  techniques.  We  have  compared  [155]  the  compression  ratio  of 
the  Lempel-Ziv  algorithms  with  the  empirical  entropy  of  the  input  string.  We  show  that 
although  these  algorithms  are  optimal  according  to  the  generally  accepted  definition,  we  can 
find  families  of  low  entropy  strings  which  are  not  compressed  optimally.  We  then  present 
a  compression  algorithm  which  combines  Lempel-Ziv  with  Run  Length  Encoding,  and  we 
show  that  it  has  excellent  performance. 

In  [148]  we  describe  a  host  of  geometric  data  structures  and  how  they  can  be  used  to 
efficiently  store  and  query  geometric  information.  In  [124]  we  show  that  the  classic  k- D 
tree  data  structure  actually  has  efficient  polylogarithmic  worst-case  time  performance  for 
answering  approximate  nearest-neighbor  searches,  which  confirms  in  a  theoretical  analysis  a 
behavior  that  researchers  have  often  observed  in  practice. 

In  [22],  we  developed  the  first  algorithm  for  constructing  an  R-tree,  a  widely  used  data 
structure  in  practice,  with  provable  near-optimal  query  bounds.  In  [8]  we  developed  the  data 
structure  for  bulk-loading  kd-trees  and  for  performing  updates.  We  wrote  two  comprehensive 
surveys  on  range  searching  [1,  28]. 

In  [59]  we  describe  an  efficient  algorithm  for  computing  arrangements  of  curve  segments. 
In  [60]  we  solve  a  long-standing  open  problem  in  computational  geometry,  showing  a  simple 
randomized  linear-time  algorithm  for  triangulating  an  n-vertex  simple  polygon. 

3.2.3  Clustering  problems 

Clustering  a  set  of  points  into  a  few  groups  is  frequently  used  for  statistical  analysis  and 
classification  in  numerous  applications,  including  information  retrieval,  facility  location,  data 
mining,  spatial  data  bases,  data  compression,  image  processing,  astrophysics,  and  scientific 
computing. 

In  [43]  we  presented  an  time  algorithm  for  solving  the  fc-center  problem  in 

Rd,  under  any  Lp  metric.  We  also  described  a  simple  (1  +  e)-approximation  algorithm 
for  the  ^-center  problem,  with  running  time  0(n  log  A;)  +  (k/e)0(-kl  1,d\  We  also  presented 
an  n°(fcl_1/c,)-time  algorithm  for  solving  the  L-capacitated  A:-center  problem,  provided  that 
L  =  nin/k1-1^)  or  L  =  0(1). 

In  [53],  we  presented  an  0(n4/3  log5  n)-time  algorithm  for  computing  the  discrete  2- center 
of  a  set  P  of  n  points  in  the  plane;  that  is,  computing  two  congruent  disks  of  smallest  possible 
radius,  centered  at  two  points  of  P,  whose  union  covers  P. 

One  method  to  reduce  the  dimensionality  is  projective  clustering-.  Given  a  set  S  of  n 
points  in  a  d-dimensional  space  and  two  integers  k  <  n  and  q  <  d,  find  k  (/-dimensional 
flats  hi, . . . ,  hk  and  partition  S  into  k  subsets  Si,  ...,Sk  so  that  maxi<j<fc  maxp€5,  d(p,  hi) 
is  minimized.  In  [44],  we  consider  the  following  two  instances  of  the  projective  clustering 
problem:  Given  a  set  S  of  n  points  in  and  an  integer  k  >  0,  cover  S  by  k  hyper-strips  (resp. 
hyper-cylinders)  so  that  the  maximum  width  of  a  hyper-strip  (resp.,  the  maximum  diameter 
of  a  hyper-cylinder)  is  minimized.  We  presented  efficient  approximation  algorithms  in  two 
and  three  dimensions.  In  2D  the  running  time  is  almost  linear  as  a  function  of  the  number  of 
points.  Recently,  we  developed  efficient  e-approximation  algorithms  for  projective  clustering 
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in  2D.  In  [57]  we  developed  and  implemented  practical  algorithms  for  projective  clustering 
and  applied  them  to  indexing  images. 

The  papers  [45,  46]  survey  known  geometric  algorithms  for  clustering  and  other  geometric 
optimization  problems. 

3.2.4  Data  Organization  and  Information  Visualization 

The  summarization,  storage,  and  visualization  of  large  amounts  of  data  plays  a  major  role 
in  knowledge  mining  and  decision  support  systems.  We  are  investigating  techniques  for  the 
visualization  of  relational  and  quantitative  information,  and  of  their  interplay.  Target  appli¬ 
cations  include  Web  search  engines,  command  and  control,  and  programming  environments. 
In  [126]  we  describe  a  new  geometric  data  structure  for  clustering  a  set  of  points  so  as  to 
facilitate  fast  data  mining  on  that  set.  The  data  structure,  which  we  call  the  balanced  as¬ 
pect  ratio  (BAR)  tree,  combines  the  best  features  of  octrees  and  k-d  trees,  in  that  it  achieves 
0(log  n)  height,  has  0(n)  space,  and  defines  cluster  regions  with  bounded  aspect  ratio.  Such 
a  structure  can  be  used,  for  example,  to  answer  approximate  nearest-neighbor  queries  and 
approximate  range  searching  queries  in  O(logn)  time,  and  our  experimental  analysis  shows 
that  they  run  fast  in  practice  as  well. 

In  [152]  we  show  how  to  maintain  a  planar  subdivision  (such  as  defined  by  a  Voronoi 
diagram)  so  as  to  quickly  support  point  location  queries  as  well  as  dynamic  updates  to  the 
subdivision.  Our  data  structure  achieved  O(logn)  time  updates  and  0(log2n)  time  queries. 
A  general  technique  for  dynamizing  data  structures  used  in  geometric  searching  and  graph 
drawing  is  presented  in  [115].  Optimal  data  structures  for  shortest  path  and  minimum-link 
path  queries  between  two  convex  polygons  inside  a  simple  polygonal  obstacle  are  presented 
in  [112], 

Current  work  on  information  visualization  has  focused  on  algorithms  and  systems  for 
graph  drawing.  In  [127]  we  present  a  novel  approach  for  cluster-based  drawing  of  large 
planar  graphs  that  maintains  planarity.  Our  technique  works  for  arbitrary  planar  graphs  and 
produces  a  clustering  which  satisfies  the  conditions  for  compound-planarity  (c-planarity). 
Using  the  clustering,  we  obtain  a  representation  of  the  graph  as  a  collection  of  O(logn) 
layers,  where  each  succeeding  layer  represents  the  graph  in  an  increasing  level  of  detail.  At 
the  same  time,  the  difference  between  two  graphs  on  neighboring  layers  of  the  hierarchy  is 
small,  thus  preserving  the  viewer’s  mental  map.  The  overall  running  time  of  the  algorithm 
is  0{n  log  n),  where  n  is  the  number  of  vertices  of  graph  G. 

In  [110]  we  address  the  problem  of  drawing  planar  graphs  with  circular  arcs  while  main¬ 
taining  good  angular  resolution  and  small  drawing  area.  We  present  a  lower  bound  on  the 
area  of  drawings  in  which  edges  are  drawn  using  exactly  one  circular  arc.  We  also  give 
an  algorithm  for  drawing  n- vertex  planar  graphs  such  that  the  edges  are  sequences  of  two 
continuous  circular  arcs.  The  algorithm  runs  in  0(n)  time  and  embeds  the  graph  on  the 
0(n)  x  0(n)  grid,  while  maintaining  Q(l/d(v))  angular  resolution,  where  d(y)  is  the  degree 
of  vertex  v.  Since  in  this  case  we  use  circular  arcs  of  infinite  radius,  this  is  also  the  first 
algorithm  to  simultaneously  achieve  good  angular  resolution,  small  area  and  at  most  one 
bend  per  edge  using  straight-line  segments.  Finally,  we  show  how  to  create  drawings  in 
which  edges  are  smooth  (^-continuous  curves,  represented  by  a  sequence  of  at  most  three 
circular  arcs. 
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In  [128]  we  develop  an  efficient  algorithm  for  displaying  large  graphs.  In  [131]  we  give 
a  fast  method  for  drawing  large  graphs  using  force-directed  methods.  In  [153]  we  describe 
a  linear-time  algorithm  for  drawing  n-node  planar  graphs  in  an  0(n)  x  0(n)  grid  with  few 
bends  per  edge. 

In  [106]  we  describe  a  fast  data  structure  for  performing  range  queries  in  tree  cross 
products.  In  [142]  we  describe  a  new  data  structure  for  ordered  dictionaries,  which  we  call 
the  energy-balanced  trees.  These  trees  achieve  logarithmic  amortized  update  times  and 
logarithmic  worst-case  search  time,  but  use  only  partial  rebuilding  to  achieve  balance;  they 
do  not  use  rotations. 

3.3  Methodologies 

As  mentioned  above,  earlier  research  on  geometric  computing  referred  to  a  convenient  compu¬ 
tation  model,  borrowed  from  traditional  algorithmic  research  and  based  on  the  assumptions 
of  unbounded  internal  memory  and  exact  real-number  arithmetics.  This  section  described 
our  efforts  at  correcting  the  arising  shortcomings. 

3.3.1  Algorithmic  Robustness 

Robustness  has  been  a  central  issue  throughout  our  project. 

In  [159]  we  have  elaborated  on  the  exact  computation  paradigm  and  formalized  the  notion 
of  degree  of  a  geometric  algorithm,  as  a  worst-case  quantification  of  the  precision  (number  of 
bits)  to  which  arithmetic  calculation  have  to  be  executed  in  order  to  guarantee  topological 
correctness.  Such  criterion  characterizes  the  computational  cost  of  topological  correctness 
of  the  output  for  the  case  of  noiseless  input  data,  and  therefore  spells  out  the  degree  as  an 
important  parameter  of  algorithmic  design.  Of  course,  concrete  situations  are  considerably 
more  complicated.  First  of  all,  input  data  may  be  noisy,  i.e.,  affected  by  errors  inherent 
in  physical  measurements  or  expressing  the  limited  accuracy  of  previous  computations.  In 
such  cases,  the  topology  consistent  with  the  data  is,  in  general,  not  unique  (depending  upon 
the  magnitude  of  the  noise),  and  the  algorithm  must  produce  its  own  certificate  (multiple 
topologies,  or  a  topology  consistent  with  specified  nominal  values).  Second,  and  certainly 
not  less  important,  while  the  algorithmic  degree  specifies  the  worst-case  arithmetic  demands, 
a  substantial  majority  of  the  computations  can  be  confidently  executed  on  a  floating-point 
platform,  acting  on  floating-point  approximations  of  the  underlying  integer  operands.  Here 
again,  extreme  care  must  be  exercised  in  assuring  the  validity  of  the  approximation.  This 
task  is  known  as  the  design  of  arithmetic  filters,  which  we  are  actively  investigating,  following 
our  initial  evaluation  of  some  important  geometric  tests  [118,  119].  A  main  direction  is  the 
development  of  tight  interval  analysis  for  the  very  concrete  case  of  noisy  input  data. 

In  the  degree-driven  algorithmic  design  area,  we  have  successfully  redesigned  algorithm 
for  robust  proximity  queries  [159]  and  have  revisited  [97]  the  plane-sweep  algorithms  for 
intersecting  segments.  This  is  a  fundamental  application,  but  the  available  algorithms  are 
notorious  for  not  being  robust.  Our  analysis  shows  that  the  reported  unrobustness  is  at¬ 
tributable  to  the  deployment  of  arithmetics  inadequate  for  the  specified  calculations,  and  we 
present  algorithms  that  either  match  or  closely  approach  the  lower  bounds  on  the  arithmetic 
and  are  eminently  practical,  with  no  sacrifice  of  efficiency. 
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In  [143],  we  provide  an  improved  algorithm  for  performing  snap  rounding  of  an  arrange¬ 
ment  of  line  segments  in  the  plane,  whose  techniques  conform  with  the  criterion  of  low-degree 
design,  outlined  above. 

3.3.2  External  memory  algorithms 

The  data  sets  for  many  of  today’s  computer  applications  are  too  large  to  fit  within  the 
computer’s  internal  memory  and  must  instead  be  stored  on  external  storage  devices  such  as 
disks.  A  major  performance  bottleneck  can  be  the  input/output  communication  (or  I/O) 
between  the  external  and  internal  memories.  We  have  written  the  first  survey  articles  on  the 
field  [179,  180,  139],  in  which  we  give  a  detailed  overview  of  the  state  of  the  art. 

Much  of  our  work  on  external  memory  algorithms  has  concentrated  on  geometric  problems 
with  immediate  applications  in  spatial  databases  and  geographic  information  systems  (GIS). 
Modern  GIS  systems,  and  especially  military  systems,  are  often  very  data-intensive  and  thus 
they  require  good  use  of  external  memory  techniques.  We  discuss  GIS  systems  further  in 
Section  3.5.1. 

During  the  year  we  have  intensified  our  work  on  implementing  our  external  memory  algo¬ 
rithms  in  the  TPIE  system  [177].  TPIE  (Transparent  Parallel  I/O  Environment)  is  designed 
to  allow  programmers  to  write  applications  that  make  efficient  use  of  I/O  resources.  In  [65]  we 
thus  not  only  developed  a  powerful  theoretical  framework  for  designing  efficient  algorithms 
for  large-scale  batched  dynamic  problems,  but  we  also  implemented  one  of  our  algorithms 
in  TPIE  and  showed  that  for  large  problem  sizes  it  greatly  outperforms  other  well-known 
algorithms.  We  have  also  continued  this  work  and  showed  that  our  TPIE  implementation 
of  a  theoretically  developed  algorithm  for  the  extremely  important  spatial  database  prob¬ 
lem  called  spatial  join  performs  as  well  as  the  state-of-the-art  algorithm  developed  in  the 
spatial  database  community.  Unlike  previously  known  implementations  our  implementation 
also  performs  well  on  skewed  or  higher  dimensional  data.  We  have  started  our  previously 
proposed  extension  of  TPIE  to  include  commonly  used  external  data  structures  and  are  very 
encouraged  by  the  results  obtained  so  far.  We  have  implemented  a  version  of  the  so  called 
R-trees — the  data  structure  which  has  emerged  as  an  efficient  and  effective  indexing  method 
for  spatial  data — and  demonstrated  how  the  time  spent  on  construction  such  structures  can 
be  dramatically  improved  using  external  memory  techniques  [63] . 

Many  data  sets  to  be  sorted  consist  of  a  limited  number  of  distinct  keys.  Sorting  such  data 
sets  can  be  thought  of  as  bundling  together  identical  keys  and  having  the  bundles  placed  in 
order;  we  therefore  denote  this  as  bundle  sorting.  In  [161],  we  described  an  efficient  algorithm 
for  bundle  sorting  in  external  memory.  We  show  that  our  algorithm  is  optimal  by  proving 
a  matching  lower  bound  for  bundle  sorting.  The  improved  running  time  of  bundle  sorting 
over  general  sorting  can  be  significant  in  practice,  as  demonstrated  by  experimentation. 
An  important  feature  of  the  new  algorithm  is  that  it  is  executed  “in-place” ,  requiring  no 
additional  disk  space. 

Computing  multidimensional  aggregates  in  high  dimensions  is  a  performance  bottleneck 
for  many  OLAP  applications.  Obtaining  the  exact  answer  to  an  aggregation  query  can  be 
prohibitively  expensive  in  terms  of  time  and/or  storage  space  in  a  data  warehouse  environ¬ 
ment.  It  is  advantageous  to  have  fast,  approximate  answers  to  OLAP  aggregation  queries. 
In  [181],  we  presented  a  novel  method  based  on  wavelets  that  provides  approximate  answers 
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to  high-dimensional  OLAP  aggregation  queries  in  massive  sparse  data  sets  in  a  time-efficient 
and  space-efficient  manner.  In  [162],  we  consider  the  case  in  which  the  data  are  updated 
dynamically  on  an  ongoing  basis.  Experiments  on  real  data  show  that  our  method  provides 
significantly  more  accurate  results  for  typical  OLAP  aggregation  queries  than  other  efficient 
approximation  techniques  such  as  random  sampling. 

In  [66],  we  settled  several  longstanding  open  problems  in  theory  of  indexability  and 
external  orthogonal  range  searching.  In  the  first  part  of  the  paper,  we  apply  the  theory  of 
indexability  to  the  problem  of  two-dimensional  range  searching.  We  show  that  the  special 
case  of  3-sided  querying  can  be  solved  with  constant  redundancy  and  access  overhead.  From 
this,  we  derive  indexing  schemes  for  general  4-sided  range  queries  that  exhibit  an  optimal 
tradeoff  between  redundancy  and  access  overhead.  In  [69],  we  derived  algorithms  with  similar 
performance  for  the  special  case  of  stabbing  queries,  which  is  needed  for  dynamic  interval 
maintenance. 

In  [176],  we  present  a  new  approach  to  designing  data  structures  for  the  important 
problem  of  external-memory  range  searching  in  two  and  three  dimensions.  We  construct 
data  structures  for  answering  range  queries  in  (3 ((log  log  log#  N )  log#  N  +  K/B)  I/O  oper¬ 
ations,  where  K  is  the  output  size.  We  base  our  data  structures  on  the  novel  concept  of 
^-approximate  boundaries,  which  are  manifolds  that  partition  space  into  regions  based  on 
the  output  size  of  queries  at  points  within  the  space. 

We  consider  the  problem  of  devising  external  memory  algorithms  whose  memory  alloca¬ 
tions  can  change  dynamically  and  unpredictably  at  run-time.  In  [94],  we  presented  a  simple 
and  natural  dynamic  memory  allocation  model.  We  define  memory-adaptive  external  mem¬ 
ory  algorithms  and  specify  what  is  needed  for  them  to  be  dynamically  optimal.  Using  novel 
techniques,  we  design  and  analyze  dynamically  optimal  memory-adaptive  algorithms  for  the 
problems  of  sorting,  permuting,  FFT,  permutation  networks,  (standard)  matrix  multipli¬ 
cation  and  LU  decomposition.  The  lower  bound  proof  techniques  for  sorting  and  matrix 
multiplication  are  fundamentally  distinct  techniques,  and  they  are  invoked  by  most  other 
external  memory  lower  bounds;  hence  we  anticipate  that  the  techniques  presented  here  will 
apply  to  many  external  memory  problems. 

In  [92],  we  consider  a  cache  shared  by  several  concurrently  running  application  processes 
and  propose  a  provably  efficient  application-controlled  global  strategy  for  the  shared  cache. 
Using  future  information  implicitly  in  the  form  of  good  decisions  by  application  processes, 
we  are  able  to  break  through  the  Hk  lower  bound  on  competitive  ratio  proved  for  classical 
paging  for  a  k- sized  cache.  For  a  siz e-k  cache  shared  by  P  application  processes  that  always 
make  good  cache  replacement  decisions,  we  develop  an  online  application-controlled  paging 
algorithm  with  a  competitive  ratio  of  2Hp-\  +  2.  Typically,  P  is  much  smaller  than  k, 
perhaps  by  several  orders  of  magnitude.  Our  competitive  ratio  improves  upon  the  2P  +  2 
competitive  ratio  achieved  by  Cao  et  al.  We  show  for  this  problem  that  no  online  algorithm 
A  can  have  a  competitive  ratio  better  than  i7p_i  even  if  the  application  processes  aiding  A 
have  perfect  knowledge  of  individual  request  sequences.  Our  results  are  with  respect  to  a 
worst-case  interleaving  of  the  individual  request  sequences  of  the  P  applications.  We  also 
consider  other  realistic  notions  of  fairness. 

In  [93],  we  provide  a  competitive  analysis  framework  for  online  prefetching  and  buffer 
management  algorithms  in  parallel  I/O  systems,  using  a  read-once  model  of  block  references. 
This  has  widespread  applicability  to  key  I/O-bound  applications  such  as  external  merging 
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and  concurrent  playback  of  multiple  video  streams.  Two  realistic  lookahead  models,  global 
lookahead  and  local  lookahead,  are  defined.  Algorithms  NOM  and  GREED  based  on  these 
two  forms  of  lookahead  are  analyzed  for  shared  buffer  and  distributed  buffer  configurations, 
both  of  which  occur  frequently  in  existing  systems.  An  important  aspect  of  our  work  is 
that  we  show  how  to  implement  both  the  models  of  lookahead  in  practice  using  the  simple 
techniques  of  forecasting  and  flushing. 

In  [6],  we  show  how  to  preprocess  a  set  S  of  points  in  d-dimensional  Euclidean  space 
to  get  an  external  memory  data  structure  that  efficiently  supports  linear-constraint  queries. 
Each  query  is  in  the  form  of  a  linear  constraint  a  •  x  <  b;  the  data  structure  must  report 
all  the  points  of  S  that  satisfy  the  query.  (This  problem  is  called  Halfspace  range  searching 
in  the  computational  geometry  literature.)  Our  goal  is  to  minimize  the  number  of  disk 
blocks  required  to  store  the  data  structure  and  the  number  of  disk  accesses  (I/Os)  required 
to  answer  a  query.  For  d  =  2,  we  present  the  first  near-linear  size  data  structures  that 
can  answer  linear-constraint  queries  using  an  optimal  number  of  I/Os.  We  also  present  a 
linear-size  data  structure  that  can  answer  queries  efficiently  in  the  worst  case.  We  combine 
these  two  approaches  to  obtain  tradeoffs  between  space  and  query  time.  Finally,  we  show 
that  some  of  our  techniques  extend  to  higher  dimensions. 

3.4  Emerging  Technologies 

3.4.1  Kinetic  Data  Structures 

In  applications  like  virtual  reality  and  dynamic  simulations,  specified  input  objects  move  or 
deform  continuously.  In  this  context,  the  questions  of  computing  some  attribute  is  replaced 
by  that  of  maintaining  it  over  time.  We  studied  several  problems  involving  motion  using 
the  approach  of  kinetic  data  structures  (KDS)  developed  at  Stanford.  Most  of  the  work  has 
been  done  in  collaboration  with  the  MURI  center  at  Stanford. 

A  binary  space  partition  is  an  essential  structure  used  in  particular  in  visualization.  Its 
maintenance  over  time  as  objects  move  is  important  for  interactive  simulations.  In  [34],  we 
describe  the  first  known  algorithm  for  efficiently  maintaining  a  Binary  Space  Partition  (BSP) 
for  n  continuously  moving  segments  in  the  plane.  Under  reasonable  assumptions  on  the  mo¬ 
tion,  we  show  that  the  total  number  of  times  the  BSP  changes  is  0(n2),  and  that  we  can 
update  the  BSP  in  O(logn)  expected  time  per  change.  We  also  consider  the  problem  of  con¬ 
structing  a  BSP  for  n  triangles  in  three-dimensional  Euclidean  space.  We  present  a  random¬ 
ized  algorithm  that  constructs  a  BSP  of  expected  size  0(n2)  in  0(n2  log2  n)  expected  time. 
We  also  describe  a  deterministic  algorithm  that  constructs  a  BSP  of  size  0((n  +  k)  log  n) 
and  height  O(logn)  in  0((n  +  k)  log2  n)  time,  where  k  is  the  number  of  intersection  points 
between  the  edges  of  the  projections  of  the  triangles  onto  the  xy- plane.  We  later  extended 
this  result  to  maintaining  a  BSP  of  moving  polyhedral  objects  [29].  We  also  proved  lower 
bounds  for  kinetic  BSPs  and  triangulations  [19]. 

In  [20,  130],  we  presented  kinetic  data  structures  for  detecting  collisions  between  a  set  of 
polygons  that  are  not  only  moving  continuously  but  whose  shapes  can  also  change  continu¬ 
ously  with  time.  We  construct  a  planar  subdivision  of  the  common  exterior  of  the  polygons, 
called  pseudo-triangulation,  which  certifies  their  disjointness.  See  also  [129]. 

We  applied  the  kinetic  setting  to  a  very  different  set  of  problems  in  parametric  optimiza- 
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tion.  Several  questions  that  involve  finding  an  optimum  tradeoff  between  two  quantities 
(e.g.  cost  and  reliability)  can  be  stated  as  parametric  optimization  problem.  We  studied  the 
parametric  minimum  spanning  tree  problem  and  obtained  substantial  improvements  over 
previous  known  solutions  [27]. 

In  [33]  we  developed  a  kinetic  data  structure  for  maintaining  various  extants  of  moving 
points,  including  diameter,  width,  and  minimum  bounding  rectangle.  It  turns  out  that 
maintaining  any  of  these  extants  exactly  is  expensive,  so  we  recently  developed  considerably 
faster  algorithms  for  maintaining  them  approximately  [35].  In  [5],  we  proposed  various 
indexing  schemes  to  store  a  set  S  of  N  points  in  the  plane,  each  moving  along  a  linear 
trajectory,  so  that  a  query  of  the  following  form  can  be  answered  quickly:  Given  a  rectangle 
R  and  a  real  value  t  (called  time  stamp),  report  all  points  of  S  that  lie  inside  R  at  time  t. 
We  propose  algorithms  that  use  kinetic  data  structures  as  well  time-space  framework.  We 
obtain  a  tradeoff  between  the  query  time  and  the  time  spent  in  evolving  the  data  structure 
as  the  points  move.  We  also  present  an  indexing  scheme  whose  query  time  depends  on  the 
value  of  t,  the  time  stamp  of  the  query  t.  In  [9]  we  extended  these  techniques  to  obtain  a 
trade-off  between  query  time  and  the  time  spent  in  updating  the  structure.  In  recent  work, 
we  developed  a  novel  method  for  maintaining  an  R-tree  on  a  set  of  moving  points  in  the 
plane.  An  R-tree  is  a  commonly  used  data  structure  for  answering  various  geometric  queries. 

3.4.2  Elements  of  a  Geometric  Library  (GeomLib) 

The  GeomLib  project  addresses  the  important  objective  of  developing  an  easy  to  use,  reliable, 
open  library  of  robust  and  efficient  geometric  algorithms.  Recent  object-oriented  design 
concepts  such  as  design  patterns  and  algorithm  abstraction  are  extensively  used  throughout 
the  entire  project.  In  the  design  phase,  we  have  taken  into  account  the  experience  of  other 
similar  efforts,  such  as  the  Library  of  Efficient  Data  structures  and  Algorithms  (LEDA)  and 
the  more  recent  Computational  Geometry  Algorithms  Library  (CGAL). 

•  To  provide  researchers  in  computational  geometry  with  a  framework  for  algorithm 
engineering,  with  a  specific  emphasis  on  geometric  computing.  In  this  context,  Geom¬ 
Lib  will  be  typically  used  for  rapid  prototyping  and  experimental  studies  of  geometric 
algorithms. 

•  To  make  computational  geometry  results  available  to  the  users  in  other  areas,  such  as 
robotics,  geographic  information  systems  (see  Section  3.5.1),  mechanical  engineering, 
computer  graphics,  etc. 

The  principles,  architecture,  and  design  of  GeomLib  are  presented  in  [172].  This  paper 
also  demonstrates  the  applicability  of  algorithm  engineering  and  software  design  concepts  to 
geometric  computing  through  a  vertical  case  study  on  the  implementation  of  planar  point 
location  algorithms  within  GeomLib  (see  also  http://www.jdsl.org/src/pointloc/).  The  design 
and  implementation  of  the  core  data  structures  support  for  GeomLib  is  presented  in  [144, 
145,  165].  Related  work  on  checkers  for  verifying  the  correctness  of  geometric  structures 
such  as  convex  polytopes  and  planar  subdivisions  are  presented  in  [117].  Related  work  on 
algorithm  engineering  and  applications  to  computer  science  education  is  presented  in  [95, 
136,  151], 
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A  major  milestone  in  the  development  of  GeomLib  has  been  the  release  of  version  2.0 
the  combinatorial  component  of  GeomLib  under  the  name  of  JDSL  in  August  2000  [170] 
(see  http://www.jdsl.org/).  Six  months  after  its  release,  JDSL  2.0  has  been  downloaded  by 
more  than  2,500  users  worldwide.  Also,  more  than  100  users  have  expressed  interest  in 
collaborating  to  the  extension  of  the  library. 

3.4.3  Geometric  Internet  Computing  (GeomNet) 

In  [78]  we  present  GeomNet,  a  system  for  performing  distributed  geometric  computing  over 
the  Internet.  The  mail  goal  of  the  system  is  facilitating  software  evaluation  and  adoption, 
as  well  as  serving  for  many  other  purposes,  through  providing  easy  (Internet)  access  to  a 
variety  of  complex  geometric  algorithms.  We  describe  the  architecture  of  GeomNet  as  a  series 
of  layers  and  discuss  the  connections  between  them.  We  also  provide  several  examples  of 
geometric  algorithms  that  our  system  already  supports.  Application  domains  for  GeomNet 
include  collaborative  research,  distance  education,  and  software  marketing. 

We  have  developed  two  prototypes  of  geometric  Web  computing  services  under  GeomNet: 
the  Mocha  algorithm  animation  system  [72,  71,  70,  73],  and  the  Graph  Drawing  Server  [100], 

Mocha  supports  interactive  animations  over  the  Web  of  geometric  algorithms.  It  employs 
a  client-server  architecture  that  optimally  partitions  the  software  components  of  a  typical 
algorithm  animation  system,  and  leverages  the  power  of  the  Java  language.  Mocha  has 
high  levels  of  security,  protects  the  algorithm  code,  places  a  light  communication  load  on 
the  Internet,  and  allows  users  with  limited  computing  resources  to  access  animations  of 
computationally  expensive  algorithms.  The  user  interface  combines  fast  responsiveness  and 
user  friendliness  with  hypertext  narratives.  The  architecture  of  Mocha  has  been  successfully 
adopted  by  researchers  worldwide. 

The  Graph  Drawing  Server  can  be  used  for  drawing  graphs  from  user-applications,  study¬ 
ing  and  comparing  graph  drawing  algorithms,  translating  between  the  formats  for  describing 
graphs  and  their  drawings,  creating  a  database  of  graphs  occurring  in  user-applications,  an 
providing  demonstrations  in  an  educational  setting.  The  service  allows  the  user  either  to 
specify  the  graph  in  a  formal  way  (HTML  format)  or  to  interactively  draw  the  input  graph 
on  the  screen. 

3.5  Applications 

3.5.1  Geographic  Information  Systems 

Geographic  information  systems  (GIS),  which  consist  of  spatial  databases  for  storing  geo¬ 
metric  mapping  information,  offer  important  applications  for  our  research.  We  have  been 
collaborating  on  many  geometric  problems  that  arise  in  GIS  with  geometric,  geographic,  and 
GIS  researchers  worldwide,  including  those  at  the  Nicholas  School  of  Environment  at  Duke. 

The  potential  and  use  of  Geographic  Information  Systems  (GIS)  is  rapidly  increasing 
because  of  the  increasing  availability  of  massive  amounts  of  geospatial  data  from  projects 
like  NASA’s  Earth  Observing  System  and  Earth  Science  Enterprise.  However,  the  use  of 
these  massive  datasets  also  exposes  scalability  problems  with  existing  GIS  algorithms.  These 
scalability  problems  arise  from  the  fact  that  most  GIS  algorithms  have  been  designed  to 
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minimize  internal  computation  time,  whereas  I/O  communication  often  is  the  bottleneck 
when  processing  massive  amounts  of  data.  In  [67,  62,  61],  we  consider  I/O-efficient  algorithms 
for  problems  on  grid-based  terrains.  Detailed  grid-based  terrain  data  is  rapidly  becoming 
available  for  much  of  the  earth’s  surface.  We  demonstrate  the  practical  merits  of  our  work 
by  comparing  the  empirical  performance  of  our  new  algorithm  for  the  flow  accumulation 
problem  with  that  of  the  previously  best  known  algorithm.  Our  experiments  show  that 
while  our  new  algorithm  scales  nicely  with  dataset  size,  the  previously  known  algorithm 
“breaks  down”  once  the  size  of  the  dataset  becomes  bigger  than  the  available  main  memory. 

Most  spatial  join  algorithms  either  assume  the  existence  of  a  spatial  index  structure  that 
is  traversed  during  the  join  process,  or  solve  the  problem  by  sorting,  partitioning,  or  on-the- 
fly  index  construction.  In  [64],  we  developed  a  simple  plane-sweeping  algorithm  that  unifies 
the  index-based  and  non-index  based  approaches.  This  algorithm  processes  indexed  as  well 
as  non-indexed  inputs,  extends  naturally  to  multi-way  joins,  and  can  be  built  easily  from  a 
few  standard  operations.  We  present  the  results  of  a  comparative  study  of  the  new  algorithm 
with  several  index-based  and  non-index  based  spatial  join  algorithms.  We  consider  a  number 
of  factors,  including  the  relative  performance  of  CPU  and  disk,  the  quality  of  the  spatial 
indexes,  and  the  sizes  of  the  input  relations.  An  important  conclusion  from  our  work  is  that 
using  an  index-based  approach  whenever  indexes  are  available  does  not  always  lead  to  the 
best  execution  time,  and  hence  we  propose  the  use  of  a  simple  cost  model  to  decide  when  to 
follow  an  index-based  approach. 

We  have  studied  a  number  of  data  structure  engineering  issues  for  storing  and  querying 
GIS  data.  We  describe  an  adaptive  scheme  [147]  for  “morphing”  a  point-location  GIS  data 
structure  so  that  it  can  adapt  to  non-uniform  query  distributions  in  a  way  that  allows 
for  faster  than  logarithmic-time  response  to  some  queries.  We  also  provide  an  efficient 
algorithm  [150]  for  performing  ray  shooting  queries  in  two-dimensional  planar  subdivisions. 

We  have  developed  a  new  theoretical  optimal  algorithm  for  the  important  GIS  problem  of 
extracting  contour  line  from  terrain  data  stored  as  Triangulated  Irregular  Networks  (TIN)  [7]. 
We  have  also  developed  new  disk  space  and  query  efficient — and  sometimes  even  provably 
optimal — data  structures  for  the  problem  of  storing  a  set  of  points  in  Rd  on  disk  such  that 
linear-constraint  queries  can  be  answered  efficiently  [6].  A  linear-constraint  query  consists 
of  finding  all  points  lying  below  a  query  hyperplane. 

We  developed  an  I/O  efficient  dynamic  point-location  data  structure  that  uses  linear 
space,  supports  fast  insertion  and  deletion  operations,  and  can  answer  point-location  queries  [4] 
In  [55]  we  developed  approximation  algorithms  for  labeling  point  features  on  maps,  and  in 
[39]  we  described  several  heuristics  for  labeling  linear  features  (e.g.  rivers,  roads)  in  a  map. 

We  considered  the  problem  of  computing  the  shortest  path  between  two  points  on  a 
given  polyhedral  surface.  This  is  a  central  problem  in  numerous  areas,  including  robotics, 
geographic  information  systems,  medical  imaging,  low-altitude  flight  simulation,  and  water- 
flow  analysis.  In  most  of  these  applications  where  efficiency  is  critical,  the  input  surfaces  are 
large  but  they  are  still  an  approximation  of  the  real  surface.  Hence,  a  simple,  efficient  algo¬ 
rithm  for  computing  an  approximate  shortest  path  is  preferable  to  an  expensive  algorithm 
that  computes  an  exact  shortest  path.  For  convex  surfaces,  we  obtained  a  linear  time  algo¬ 
rithm  that  computes  a  path  that  whose  length  is  at  most  1  +  e  that  of  the  shortest  path  [37]. 
For  non-convex  surfaces,  We  obtained  the  first  algorithms  that  compute  an  approximate 
shortest  path  in  subquadratic  time  [174].  In  [36],  we  have  developed  and  implemented  a 
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robust,  efficient  algorithm  for  computing  approximate  shortest  paths  on  a  convex  polyhedral 
surface.  Although  the  analysis  and  the  current  implementation  of  the  algorithm  work  only 
for  convex  polytopes,  many  of  ideas  extend  to  arbitrary  polyhedral  surfaces. 

In  [68],  we  develop  efficient  new  external-memory  algorithms  for  a  number  of  important 
problems  involving  line  segments  in  the  plane,  including  trapezoid  decomposition,  batched 
planar  point  location,  triangulation,  red-blue  line  segment  intersection  reporting,  and  general 
line  segment  intersection  reporting.  In  GIS  systems,  the  first  three  problems  are  useful  for 
rendering  and  modeling,  and  the  latter  two  are  frequently  used  for  overlaying  maps  and 
extracting  information  from  them. 

3.5.2  Robotics 

We  consider  the  problem  of  computing  the  shortest  path  between  two  points  on  a  given  poly¬ 
hedral  surface.  This  is  a  central  problem  in  numerous  areas,  including  robotics,  geographic 
information  systems,  medical  imaging,  low-altitude  flight  simulation,  and  water-flow  analy¬ 
sis.  In  most  of  these  applications,  a  simple,  efficient  algorithm  for  computing  an  approximate 
shortest  path  is  preferable  to  an  expensive  algorithm  that  computes  an  exact  shortest  path, 
since  the  input  surfaces  are  large,  efficiency  is  critical,  and  the  polyhedral  surface  is  typically 
an  approximation  of  the  real  surface.  All  previous  approaches  take  superquadratic  time  to 
compute  even  an  approximate  shortest  path.  We  present  the  first  algorithms  that  compute 
an  approximate  shortest  path  in  subquadratic  time.  An  undergraduate  student  is  currently 
implementing  our  algorithm. 

When  a  robot  is  too  big  to  be  assimilated  to  a  point,  the  next  approximation  is  to  consider 
the  smallest  ball  containing  it.  If  we  wish  to  obtain  a  path  for  the  ball  that  doesn’t  collide 
with  any  obstacle  in  space,  a  standard  approach  is  to  reduce  this  problem  to  that  of  planning 
the  motion  of  a  point  in  a  space  where  each  obstacle  is  grown  by  the  radius  of  the  ball.  We 
show  that  this  space  has  a  complexity  roughly  quadratic  in  the  number  of  objects  [50,  47]. 
Only  the  trivial  cubic  bound  was  known  previously. 

We  have  also  been  studying  the  optimal  motion-planning  problems  under  nonholonomic 
constraints.  We  have  developed  a  simple,  efficient  algorithm  for  computing  a  curvature- 
constrained  approximate  shortest  path  for  a  point  robot  moving  amid  planar  obstacles  [182]. 
We  have  also  developed  a  simple,  efficient  algorithm  for  computing  a  curvature-constrained 
approximate  shortest  path  for  a  point  robot  moving  inside  a  convex  polygon.  We  proved 
several  interesting  properties  of  shortest  paths,  which  we  believe  are  useful  for  navigation 
amid  planar  obstacles  [21].  In  collaboration  with  the  robotics  group  at  Stanford,  we  studied 
the  problem  of  moving  an  object  by  pushing  it  [40].  See  [2,  3,  16,  23,  30,  12,  32]  for  other 
robotics  problems,  including  assembly  design,  penetration  depth,  and  polygon  placement, 
which  we  studied. 

3.5.3  Information  Visualization 

The  visualization  of  large  amounts  of  data  plays  a  major  role  in  information  retrieval,  knowl¬ 
edge  mining  and  decision  support  systems,  which  are  critical  applications  especially  in  re¬ 
lation  to  the  digital  battlefield.  In  this  context  we  have  devoted  particular  attention  to 
geometrically  visualizing  graph  and  networks  and  we  have  thus  focused  on  graph  drawing 
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algorithms  and  systems. 

The  way  in  which  space  is  partitioned  relative  to  a  geometric  model  can  also  have  a  big 
impact  on  rendering  speed.  Thus,  as  another  step  towards  developing  efficient  rendering 
algorithms,  we  have  been  studying  binary  space  partitions  (BSP),  which  is  a  data  structure 
of  fundamental  importance  in  visualization,  with  applications  in  hidden-surface  removal,  ray 
tracing,  global  illumination,  and  shadow  generation.  All  the  lower-bound  constructions  for 
large-size  BSPs  occur  very  rarely  in  practice.  Continuing  our  work  on  constructing  BSPs 
for  orthogonal  rectangles  [31,  41],  we  studied  the  problem  of  constructing  BSPs  for  triangles 
[34],  Our  algorithms  perform  very  efficiently  particularly  when  the  triangles  model  terrain¬ 
like  surfaces.  With  the  emergence  of  VRML,  many  geometric  data  sets  contain  moving 
objects.  Motivated  by  this  fact,  we  have  developed  a  fast  algorithm  for  maintaining  the  BSP 
for  continuously  moving  triangles,  which  we  describe  in  Section  3.4.1. 

In  [98,  99],  we  address  the  visualization  of  large-scale  transportation  and  communication 
networks  and  present  a  new  edge  routing  method  that  uses  cubic  curves  in  the  plane  and  on 
the  sphere.  We  also  provide  novel  interactive  visualizations  of  air  and  train  connections  and 
of  multicast  Internet  traffic.  Experiments  on  real-world  data  sets  indicate  that  our  method 
is  computationally  fast  and  can  effectively  visualize  networks  with  7,000  nodes  and  10,000 
edges. 

In  [102],  we  introduce  a  framework  for  defining  and  validating  metrics  to  measure  the 
difference  between  two  drawings  of  the  same  graph,  and  give  a  preliminary  experimental 
analysis  of  several  simple  metrics.  In  [103],  we  formally  define  several  intuitive  ideas  of 
similarity  and  present  the  results  of  a  user  study  designed  to  evaluate  how  well  these  measures 
reflect  human  perception  of  similarity. 

In  [128]  we  develop  an  efficient  algorithm  for  displaying  large  graphs  using  geometric 
space-partitioning  techniques.  In  [131]  we  give  a  fast  method  for  drawing  large  graphs  using 
force-directed  methods. 

A  systematic  study  of  planar  drawings  has  been  conducted  with  the  goal  of  optimizing  the 
embedding  of  graphs  in  the  plane.  An  optimal  compaction  method  for  orthogonal  represen¬ 
tations  is  given  in  [104].  Convex  planar  grid  drawings  are  studied  in  [122,  123].  Planar  tree 
drawings  are  studied  in  [109,  132,  154].  A  parallel  orthogonal  drawing  technique  is  presented 
in  [171].  Bend  and  area  minimization  in  curvilinear  drawings  is  addressed  in  [153,  110]. 
In  [163],  we  give  a  linear  time  algorithm  for  computing  a  minimum-depth  embedding  of 
a  planar  graph.  A  linear-time  algorithm  for  testing  the  upward  planarity  of  single-source 
digraphs  is  presented  in  [96].  In  [134],  we  show  that  upward  planarity  testing  and  rectilin¬ 
ear  planarity  testing  are  NP-complete  problems.  In  [127]  we  present  a  novel  approach  for 
cluster-based  drawing  of  large  planar  graphs  that  maintains  planarity. 

A  system  for  the  interactive  construction  of  orthogonal  drawings  is  presented  in  [105].  An 
object-oriented  design  and  a  software  prototype  of  a  fundamental  orthogonal  graph  drawing 
technique  is  presented  in  [137] 

Three-dimensional  drawings  of  graphs  are  studied  in  [113,  116,  133,  135]. 

In  [158, 13,  121,  160],  we  explore  various  types  of  proximity  drawings,  characterize  classes 
of  graphs  that  admit,  such  drawings,  and  investigate  the  time  complexity  of  constructing 
them. 

In  [101],  we  present  an  application  of  our  graph  drawing  and  Internet  computing  tech¬ 
niques  to  computer  science  education. 
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The  book  [120]  is  recognized  as  the  authoritative  reference  in  the  field  of  graph  draw¬ 
ing  and  can  be  used  as  a  textbook  in  graduate  courses.  A  survey  of  research  on  graph 
drawing  is  given  in  [168].  Recent  developments  in  graph  drawing  research  are  overviewed 
in  [166]  and  [167],  In  [178],  we  survey  experimental  studies  on  graph  drawing  algorithms. 
Techniques  for  supporting  constraints  in  a  wide  variety  of  graph  drawing  algorithms  are 
discussed  in  [164]. 

3.5.4  Geometric  Modeling  and  Computer  Graphics 

Computer-aided  geometric  modeling  and  computer  graphics  are  important  topics  to  the 
Army  because  of  the  central  role  it  plays  in  vehicle  design,  testing,  and  battle  visualization 
and  simulation.  Of  particular  significance  is  the  role  that  BRL-CAD  plays  in  establishing 
the  vulnerability  of  various  military  assets.  We  have  therefore  carried  out  several  activities 
in  these  areas. 

In  [75]  we  describe  a  method  for  using  geometric  hashing  to  efficiently  fix  a  CAD  object 
whose  boundary  is  broken  due  to  imprecision,  erroneous  approximation  of  smooth  surfaces, 
etc.  The  algorithm  preserves  global  orientability  of  the  repaired  object.  We  continue  and 
extend  work  on  CAD  repair  in  work  that  has  been  of  significant  interest,  particularly  to 
our  colleagues  at  ARL-Aberdeen.  Specifically,  we  have  designed  a  general  framework  and 
new  model-repair  heuristic  algorithm  that  couples  polygon  merging  with  a  visualization 
system  for  repairing  CAD  models  [90].  In  addition,  we  provide  a  “bridge”  [89]  between 
CAD  systems  and  layered  manufacturing  (rapid  prototyping)  systems.  The  system  allows 
the  user  to  accommodate  CAD  files,  repair  them,  and  consolidate  multiple  objects  for  one 
machine-running  job. 

We  have  also  investigated  a  number  of  theoretical  and  applied  issues  regarding  the  general 
problem  of  storing  and  representing  sets  of  polygons  in  three-dimensional  space.  We  designed 
and  implemented  a  software  system  [74]  that  supports  storage  of  polyhedral  data  and  of 
implementing  algorithms  on  polyhedra.  A  three-dimensional  polygon  is  a  closed  polygonal 
chain  is  3D.  We  show  [81]  that  determining  its  triangulability  is  NP-complete,  and  we  show 
some  sufficient  and  necessary  conditions  for  triangulability. 

In  [31,  42,  41],  we  consider  the  practical  problem  of  constructing  binary  space  partitions 
(BSPs)  for  a  set  S  of  n  orthogonal,  non-intersecting,  two-dimensional  rectangles  in  three- 
dimensional  Euclidean  space  such  that  the  aspect  ratio  of  each  rectangle  in  S  is  at  most 
a,  for  some  constant  a  >  1.  We_present  an  0(n2V^)-time  algorithm  to  build  a  binary 
space  partition  of  size  0(n21Ae")  for  g  We  also  show  tjiat  m  0f  the  n  rectangles  in 
S  have  aspect  ratios  greater  than  a,  we  can  construct  a  BSP  of  size  for 

S  in  0(riy/m2vlosn)  time.  The  constants  of  proportionality  in  the  big-oh  terms  are  linear 
in  loga.  We  extend  these  results  to  cases  in  which  the  input  contains  non-orthogonal  or 
intersecting  objects.  We  also  developed  several  approximation  algorithms  for  surface,  curve, 
and  map  simplification,  a  problem  central  to  multiresolution  modeling  [54,  25,  114,  56]. 
We  [38]  developed  a  novel  occlusion-culling  algorithm  that  exploits  the  graphics  hardware 
and  the  geometry  of  the  input. 

In  [157],  we  explore  a  novel  and  practical  method  for  back-face  culling,  which  allows 
rendered  scenes  to  be  displayed  faster  in  a  visually-accurate  manner.  We  continued  this 
approach  in  [85]  and  [86],  where  we  explore  how  to  exploit  coherence  in  interactive  and  fly- 
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through  visual  simulations  of  a  polyhedral  model  to  display  silhouettes  in  highlight.  This 
approach  allows  us  to  focus  on  the  main  factors  that  define  the  “shape”  of  an  object  in  a 
simulation.  Specifically,  we  present  an  efficient  algorithm  for  computing  a  view-dependent 
silhouette  of  a  polyhedral  model.  We  compute  perspective-accurate  silhouettes  by  reducing 
the  problem  to  point  location  queries,  and  demonstrate  the  usefulness  of  a  fast  silhouette- 
finder  in  rendering  relatively  large  models,  simplification  and  medical  data  registration. 

In  addition,  we  have  developed  a  course  [58]  at  SIGGRAPH  that  explores  these  topics  of 
efficiency  in  the  context  of  performing  interactive  walk-throughs  in  large  geometric  databases, 
particularly  those  of  interest  in  defense  settings. 

3.5.5  Geometric  Pattern  Matching  and  Computational  Metrology 

Determining  if  a  set  of  points  (taken  as  samples  from  some  manufactured  part)  conforms  to 
a  specification  is  essentially  a  geometric  pattern  matching  problem.  Defining  and  analyzing 
geometric  models  is  often  a  first  step  in  the  manufacturing  process,  with  the  final  step  being 
that  of  matching  that  model  to  a  set  of  sensor  data  (for  either  computational  metrology  or 
for  automatic  target  recognition). 

We  initiated  pioneering  work  that  precisely  defined  the  point-set  pattern  matching  prob¬ 
lem  from  an  object-space  point  of  view,  and  we  gave  several  efficient  deterministic  algorithms 
for  solving  geometric  pattern  matching  problems  in  this  framework  [111].  Moreover,  we 
provide  several  efficient  exact  and  approximation  algorithms  [125]  for  solving  object-space 
geometric  pattern  matching  problems  that  arise  in  computational  metrology  applications. 

In  addition,  we  define  a  new  polygon-offset  distance  function  and  show  how  to  efficiently 
construct  Voronoi  diagrams  using  this  measure  [83].  This  distance-function  is  motivated  by 
geometric  pattern  matching  problems  involving  polygonal  shapes,  and  we  provide  a  number 
of  applications  [80]  of  our  polygon-offset  distance  function  definition  and  Voronoi  diagrams 
defined  by  this  distance  to  problems  in  computational  metrology  and  robot  motion  plan¬ 
ning.  We  produced  a  video  [79]  published  in  the  video  proceedings  of  the  Symposium  on 
Computational  Geometry  that  was  an  animation  of  our  algorithm  for  fitting  a  polygon  to 
a  set  of  points,  as  in  a  computational  metrology  application,  by  using  a  new  “distance” 
measure  that  is  defined  by  computing  offsets  from  a  reference  polygon.  We  also  describe  two 
improved  algorithms  [84]  for  translating  a  convex  polygon  to  contain  a  maximum  number 
of  points,  which  has  applications  in  computer  vision  and  geometric  optimization,  as  well  as 
provides  a  robust  algorithm  [91]  for  partial  matching  between  surfaces  and  volumes,  where 
the  existence  of  no  predefined  features  is  assumed. 

We  present  in  [87]  an  efficient  0(n  +  l/e4-5)-time  algorithm  for  computing  a  (1  +  e)- 
approximation  of  the  minimum- volume  bounding  box  of  n  points  in  5?3.  We  also  present  a 
simpler  algorithm  (for  the  same  purpose)  whose  running  time  is  0(nlogn  +  n/e3). 

In  [11,  77,  80]  we  give  solutions  to  some  constrained  annulus  placement  problems  for  offset 
polygons  and  scaled  polygons.  Such  matching  problems  show  how  closely  a  set  of  sample 
points  matches  a  circular  or  polygonal  boundary.  The  goal  is  to  find  the  smallest  annulus 
region  of  a  circle  or  polygon  containing  a  set  of  points  subject  to  the  following  constraints: 
we  fix  the  inner  (resp.,  outer)  polygon  of  the  annulus  and  minimize  the  annulus  region  by 
minimizing  the  outer  offset  (resp.,  maximizing  the  inner  offset).  We  also  solve  a  problem 
that  can  be  viewed  as  a  special  case  of  the  first  problem:  finding  the  smallest  translated 
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copy  of  a  polygon  containing  an  entire  point  set.  For  all  of  these  problems,  we  solve  for  the 
cases  where  smallest  and  largest  are  defined  by  either  offsetting  the  polygon  or  scaling  it. 
In  [14]  we  developed  algorithms  for  computing  the  smallest  cylinder  enclosing  a  point  set, 
and  in  [17]  we  developed  an  approximation  algorithm  for  fitting  a  cylinder  through  a  set  of 
points. 

The  3sum  problem  represents  a  class  of  problems  conjectured  to  require  fl(n2)  time  to 
solve,  where  n  is  the  size  of  the  input.  Given  two  polygons  P  and  Q  in  the  plane,  we  show 
in  [88]  that  some  variants  of  the  decision  problem,  whether  there  exists  a  transformation  of 
P  that  makes  it  contained  in  Q,  are  3sum-hard.  We  also  show  that  the  geometric  pattern 
matching  problem  of  finding  the  translation  in  the  plane  that  minimizes  the  HausdorfF  dis¬ 
tance  between  two  segment  sets  is  3sum-hard.  Nevertheless,  for  two  sets  of  points  P  and 
Q,  we  have  designed  several  practical  algorithms  [146]  for  approximately  matching  P  to  a 
subset  of  Q  under  the  HausdorfF  distance  subject  to  rotations,  translations,  and  combined 
rotations  and  translations. 

In  [51]  we  developed  efficient  algorithms  for  exact  pattern  matching  allowing  translation 
and  rotation.  The  analysis  of  our  algorithm  relies  on  a  novel  technique  that  bounds  the 
number  of  simplices  spanned  by  a  set  of  points  that  are  congruent  to  a  given  simplex. 
We  also  developed  exact  and  approximation  algorithms  for  minimum-weight  matching,  a 
fundamental  problem  in  shape  matching  [173,  175]. 

In  the  same  line  of  research,  and  in  close  collaboration  with  a  foremost  industry  in  the 
field  (Marh/Federal  of  Providence,  RI),  we  have  developed  an  entirely  new  and  eminently 
practical  approach  to  a  central  problem  in  metrology,  the  determination  of  the  minimum 
zone  cylinder  [119].  The  minimum  zone  cylinder  of  a  set  of  points  in  three  dimensions  is 
the  cylindric  crown  defined  by  a  pair  of  coaxial  cylinders  with  minimal  radial  separation 
(width).  In  the  context  of  tolerancing  metrology,  the  set  of  points  is  nominally  cylindrical, 
i.e.,  the  points  are  known  to  lie  in  close  proximity  of  a  known  reference  cylinder.  Using 
approximations  which  are  valid  only  in  the  neighborhood  of  the  reference  cylinder,  we  can  get 
a  very  good  approximation  of  the  minimum  zone  cylinder.  The  process  provides  successive 
approximations,  and  each  iteration  involves  the  solution  of  a  linear  programming  problem 
in  six  dimensions.  The  error  between  the  approximation  and  the  optimal  solution  converges 
very  rapidly  (typically  in  three  iterations  in  practice)  down  to  a  limit  error  of  ^  (  where 
loq  is  the  width  and  R  is  the  external  radius  of  the  zone  cylinder). 

3.5.6  Layered  Manufacturing 

An  important  problem  in  layered  manufacturing  is  the  choice  of  a  good  build  direction,  which 
influences  the  quality  and  the  cost  of  manufacturing  the  object.  In  [26],  we  present  efficient 
algorithms  for  computing  a  build  direction  that  optimizes  the  total  area  of  faces  that  are  in 
contact  with  supports. 

4  Technology  Transfer 

Our  research  project  has  motivated  the  preliminary  establishment  of  various  collaborations 
worldwide.  Also,  attention  has  been  already  given  by  industry  to  the  potential  commercial 
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applications  of  our  work.  The  collaborations  will  provide  “real”  users  of  our  system  proto¬ 
types  and  will  leverage  various  research  efforts  funded  by  US  and  European  agencies.  The 
outreach  of  our  Center  for  Geometric  Computing  is  directed  to  both  industries  and  U.S. 
Army  (or  Department  of  Defense)  organizations. 

With  regard  to  industry,  we  have  engaged  in  collaborative  research  with  Marh/Federal 
in  the  area  of  computational  metrology  and  with  AT&T  Research  Laboratories  on  geometric 
techniques  for  storing  and  displaying  large  calling  network  graphs. 

Our  outreach  to  U.S.  Army  research  projects  has  primarily  been  through  the  Army 
Research  Lab  in  Aberdeen,  although  we  have  also  been  meeting  and  discussing  research 
problems  with  people  at  ARL-Adelphi.  We  have  also  conducted  discussions  with  scientists 
with  other  Department  of  Defense  organizations,  particularly  DARPA  and  NSA.  Specifically, 
our  discussions  have  been  directed  at  studying  the  feasibility  of  using  our  geometric  data 
structures  library  for  Internet  security  infrastructure  tasks,  such  as  certificate  repository,  as 
well  as  mechanisms  for  adding  security  to  geometric  collaboration  tasks,  such  as  exist  in  our 
GeomNet  system. 

Our  efforts  are  also  having  an  impact  internationally.  For  example,  our  project  on  soft¬ 
ware  development  and  Internet  computing  is  likely  to  have  an  impact  on  the  development 
of  the  CGAL  library  of  geometric  algorithms,  a  major  European  research  project  with  seven 
sites.  International  collaborations  with  the  INRIA  French  site  (Jean-Daniel  Boissonnat,  di¬ 
rector),  the  Saarbriicken  German  site  (Kurt  Mehlhorn,  director),  and  the  Israeli  site  (Micha 
Sharir,  director)  have  been  activated  and  materialized  in  productive  visits. 


